Approximating the Matrix Exponential of an Advection-Diffusion Operator Using the Incomplete Orthogonalization Method
نویسنده
چکیده
In this paper we give first results for the approximation of eb, i.e. the matrix exponential times a vector, using the incomplete orthogonalization method. The benefits compared to the Arnoldi iteration are clear: shorter orthogonalization lengths make the algorithm faster and a large memory saving is also possible. For the case of three term orthogonalization recursions, simple error bounds are derived using the norm and the field of values of the projected operator. In addition, an a posteriori error estimate is given which in numerical examples is shown to work well for the approximation. In the numerical examples we particularly consider the case where the operator A arises from spatial discretization of an advection-diffusion operator.
منابع مشابه
Galerkin Method for the Numerical Solution of the Advection-Diffusion Equation by Using Exponential B-splines
In this paper, the exponential B-spline functions are used for the numerical solution of the advection-diffusion equation. Two numerical examples related to pure advection in a finitely long channel and the distribution of an initial Gaussian pulse are employed to illustrate the accuracy and the efficiency of the method. Obtained results are compared with some early studies.
متن کاملPreconditioned Generalized Minimal Residual Method for Solving Fractional Advection-Diffusion Equation
Introduction Fractional differential equations (FDEs) have attracted much attention and have been widely used in the fields of finance, physics, image processing, and biology, etc. It is not always possible to find an analytical solution for such equations. The approximate solution or numerical scheme may be a good approach, particularly, the schemes in numerical linear algebra for solving ...
متن کاملThe Structure of Bhattacharyya Matrix in Natural Exponential Family and Its Role in Approximating the Variance of a Statistics
In most situations the best estimator of a function of the parameter exists, but sometimes it has a complex form and we cannot compute its variance explicitly. Therefore, a lower bound for the variance of an estimator is one of the fundamentals in the estimation theory, because it gives us an idea about the accuracy of an estimator. It is well-known in statistical inference that the Cram&eac...
متن کاملPositivity-preserving nonstandard finite difference Schemes for simulation of advection-diffusion reaction equations
Systems in which reaction terms are coupled to diffusion and advection transports arise in a wide range of chemical engineering applications, physics, biology and environmental. In these cases, the components of the unknown can denote concentrations or population sizes which represent quantities and they need to remain positive. Classical finite difference schemes may produce numerical drawback...
متن کاملSpace-time radial basis function collocation method for one-dimensional advection-diffusion problem
The parabolic partial differential equation arises in many application of technologies. In this paper, we propose an approximate method for solution of the heat and advection-diffusion equations using Laguerre-Gaussians radial basis functions (LG-RBFs). The results of numerical experiments are compared with the other radial basis functions and the results of other schemes to confirm the validit...
متن کامل